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We propose a novel method to probe primordial gravitational waves by means of primordial black 
holes (PBHs). When the amplitude of primordial tensor perturbations on comoving scales much 
smaller than those relevant to Cosmic Microwave Background is very large, it induces scalar pertur¬ 
bations due to second-order effects substantially. If the amplitude of resultant scalar perturbations 
becomes too large, then PBHs are overproduced to a level that is inconsistent with a variety of 
existing observations constraining their abundance. This leads to upper bounds on the amplitude of 
initial tensor perturbations on super-horizon scales. These upper bounds from PBHs are compared 
with other existing bounds. 

PACS numbers: 


Introduction.- Stochastic gravitational wave back¬ 
ground (SGWB) on a wide range of scales is thought 
to have been generated in the early universe. SGWB on 
the larTOst observable scales have been investigated by 
Planck Q and BICEP2 Q. SGWB on smaller scales can 
be constrained by inferring the value of iVeff, the effective 
number of extra degrees of freedom of relativistic species, 
at Big Bang Nucleosynthesis (BBN) through the current 
abundance of light elements Q , or at photon decoupling 
through the anisotropy of Cosmic Microwave Background 
(CMB) 0,0. More recently SGWB on small scales has 
been constrained by CMB spectral distortions as well 
00. SGWB can also be directly constrained by gravita¬ 
tional wave detectors (see e.g. 0). BBN and CMB have 
played a major role in constraining SGWB since they are 
applicable on a wide range of scales, noting ground-based 
laser interferometers tend to target gravitational waves 
(GWs) on a relatively limited frequency range with high 
sensitivity. However, it would be worthwhile to note that 
upper bounds obtained through need an assumption 
about the number of relativistic species in the early uni¬ 
verse, as is discussed later. In addition, in obtaining BBN 
or CMB bounds we implicitly assume that any physical 
mechanisms, both known and unknown, increase N^g, 
making TVeff larger than the standard value A^eff = 3.046 
0. However, in principle it is possible to decrease N^g 
(see e.g. [IMl)- Therefore, it would be desirable to 
have another independent cosmological method to probe 
SGWB on a wide range of scales, which does not much 
depend on the assumptions mentioned above. 

The very physical mechanism we employ to probe 
SGWB is the formation of primordial black holes (PBHs), 
black holes formed in the early universe well before the 
cosmic structure formation. There are a number of pos¬ 
sible mechanisms to create PBHs (see e.g. [l^ and refer¬ 
ences therein), but one of the simple and natural mecha¬ 
nisms is the direct collapse of radiation overdensity dur¬ 
ing the radiation-dominated era, which happens when 


the density perturbation becomes order unity at the mo¬ 
ment of the horizon crossing of the perturbation [IM3. 
(See also [13, [3 updated discussions of the formation 
condition and see also [Tol - I^ for numerical simulations 
of the formation of PBHs.) There is no conclusive ev¬ 
idence for the existence of PBHs and so upper bounds 
on the abundance of PBHs on various mass scales have 
been obtained by various kinds of observations (see e.g. 
0 and references therein). These upper bounds can 
be translated into upper bounds on the power s pec trum 
of the curvature perturbation on small scales |26l [^ . 
namely, they can be used to exclude models of the early 
universe which predict too many PBHs. (Other methods 
to constrain primordial perturbations on small scales in¬ 
clude GMB spectral distortions [^, acoustic reheating 
[23,[33|, ultracompact minihalos [SI, HI] etc.) 

In this paper we use PBHs to constrain primordial 
tensor fluctuations. Large amplitude tensor perturba¬ 
tions induce scalar perturbations {induced scalar pertur¬ 
bations) due to their second-order effects, though tensor 
and scalar perturbations are decoupled at the linear level. 
If the typical amplitude of primordial GWs is too large, 
then the typical amplitude of resultant induced scalar 
perturbations becomes also too large and there appear 
too many regions where the amplitude of density per¬ 
turbation becomes order unity at the horizon crossing to 
form PBHs. Thus we can obtain upper bounds on the 
amplitude of primordial tensor modes requiring PBHs are 
not overproduced as they enter the Hubble radius [^ . 

Due to our current ignorance of the correct model 
of inflation and the subsequent thermal history of the 
universe, new upper limits on tensor perturbations on 
small scales in themselves are worthwhile. In addition, 
it makes our new upper limits still more valuable that 
there are models of the early universe which can predict 
large tensor perturbations on small scales (33i - l4^ . If a 
model predicts large tensor perturbations on small scales, 
and even larger scalar perturbations at the same time. 
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then such a model would be more severely constrained 
by PBHs generated from first-order scalar perturbations. 
Here we consider PBH formation only from second-order 
tensor perturbations, but if scalar perturbations are also 
large, PBHs are formed more, and so our bounds on ten¬ 
sor perturbations are conservative or model-independent, 
in the sense that those bounds do not depend on the am¬ 
plitude of first-order scalar perturbations on small scales. 
Note also that some of the above previous works can (or 
possibly can) predict not only large tensor perturbations 
on small scales, but also large tensor-to-scalar ratio on 
small scales, and our method is particularly useful to con¬ 
strain these types of models. 

To calculate induced scalar perturbations one may 
solve the Einstein equations for scalar perturbations in¬ 
cluding the source terms which are second-order in ten¬ 
sor perturbations, assuming some form of the initial 
power spectrum of tensor perturbations. In general, the 
statistics of induced scalar perturbations are highly non- 
Gaussian, since they are generated by second-order ten¬ 
sor perturbations. If we assume some statistical prop¬ 
erties of initial tensor perturbations, e.g. Gaussianity, 
in principle one can determine the probability density 
function (PDF) of induced scalar perturbation, which is 
necessary to calculate the abundance of PBHs given some 
initial tensor power spectrum to obtain PBH constraints 
on primordial GWs. All the details of those analyses will 
be presented elsewhere [43| . and here we present simple 
analytic estimations, which correctly reproduce the main 
results of Q, as is discussed shortly. 

PBH upper bounds on tensor perturbation.- For sim¬ 
plicity we assume the typical amplitude of tensor fluctu¬ 
ations of only a limited range of wavelengths around \p 
(say, in the range (e“^/^Ap, e^^^Ap)) is large, namely we 
consider a peak in the power spectrum around Ap. Let 
us denote by Oh the root-mean-square amplitude (RMS) 
or the typical amplitude of tensor fluctuations when the 
relevant modes under consideration are on super-horizon 
scales. Note that roughly corresponds to the dimen¬ 
sionless power spectrum Vh of primordial GWs. Sup¬ 
pose these large amplitude tensor perturbations whose 
wavelength is around Ap reenter the horizon during the 
radiation-dominated era. It generates radiation density 
perturbations 5r whose wavelength is around Ap and 
whose RMS at around the horizon crossing is roughly 
given by 175 ^ ~ cr^, simply because 5^ is generated by 
second-order GWs. As confirmed in [i^, this relation 
indeed holds well at least for the delta-function power 
spectrum. If Sr > Sr ,th 0.4 [13 , the induced radiation 
density perturbations collapse to form PBHs. The proba¬ 
bility of the formation of PBHs has to be extremely small 
to be consistent with various observations [l^ . There¬ 
fore, in order not to overproduce PBHs, one may require 
'Ph ^ ^ <^Sr ^ 0.1(5r,th ~ 0.04. More comments will be 

made later regarding this inequality. This upper bound 
is applicable on all the scales which reenter the horizon 
during the radiation-dominated era. Furthermore, this 
bound is probably conservative on scales which reenter 
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FIG. 1: Comparison of the constraints on the amplitude of the 
delta-function power spectrum eq. ©, obtained from PBH, 
BBN, CMB, PTAs, LIGO and Virgo. For more details see 
[ 43 . The PBH constraints in this figure are obtained assum¬ 
ing Gaussianity of primordial tensor perturbations and cal¬ 
culating the PDF of induced radiation perturbations numer¬ 
ically in [ 43 . The thick and dashed blue lines correspond to 
CMB constraints for adiabatic and homogeneous initial con¬ 
ditions of GWs, respectively. Note that the constraint from 
PTA is more than two orders of magnitude tighter than those 
from PBHs at fcp ~ 4 x 10®Mpc“^. The frequency correspond¬ 
ing to each kp, fp = ckpj2n, is also shown. See the text for 
more details. 

the horizon during the (early as well as late) matter- 
dominated era, since formation of black holes by direct 
collapse of primordial perturbations is easier when the 
universe is (literally or effectively) dominated by pres¬ 
sureless dust 0 . To summarize, we have 

iP;,<0.04 (PBH), (1) 

and this is applicable from the comoving wavelength of 
~ IGpc all the way down to ~ IGpc x ^ 0.3m, 
assuming the total number of e-foldings during inflation 
is sixty ( 6 l[ . See Discussion for some comments regarding 
the constraints on smallest scales. 

Our analytic estimations can reproduce the results 
of more rigorous calculations in [44| . where the delta- 
function initial tensor power spectrum is assumed: 

Vh{k) = A^kpS{k-kp). (2) 

The PBH abundance for this power spectrum is calcu¬ 
lated and constraints on as a function of kp are ob¬ 
tained, which is shown in FiglU along with other con¬ 
straints. One can see that our analytic estimations (eqs. 
(HD, (USD and ((HD) indeed explain Fig.l well. 

Comparison with other methods.- First note that the 
total energy density of radiation without the presence of 
gravitational waves can be written as follows; 

Prad(T) = (7rV30)5*r^ (3) 
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Here 5 * is the effective number of degrees of freedom of 
relativistic species and is given by, before the electron- 
positron annihilation, [3|, l45l | 

5 * = 2 +J(4 + 2fV,), (4) 

where is the effective number of degrees of freedom 
of neutrinos = 3.046. It is convenient to denote the 
presence of GWs (or possibly dark radiation) by ANgg, 
as a correction to N^, above Q . In the following we use 
ATVgw as a contribution of GWs and derive the expres¬ 
sion of it in terms of Vh, enabling observational upper 
limits on AA^eff to be translated into upper limits on 
When GWs are present, the critical density pcrit during 
the radiation-dominated era is 


Peril — Prad “b PGW ; (5) 

where prad and pew are the energy density of radia¬ 
tion and GWs respectively. After the horizon cross¬ 
ing of GWs, their energy density starts to scale as 
oc a“^, while, denoting the effective degrees of freedom 
of relativistic species in terms of entropy at tempera¬ 
ture T by gs{T), the photon temperature evolves fol¬ 
lowing gs{T)T^a^ =const. (i.e. constant entropy) and 
so PradOc g*/{a^g‘^J^) llia^g^J^) (see e.g. [H]). Let us 
define Gqw = PGw/Perit — PGw/Prad- Then, 

^GwiT) = ^Gw{Tin), (6) 

where Tin is the temperature of radiation at the moment 
of the horizon crossing of GWs whose wavelength is ^ Xp, 
and T < Tn. At the epoch of BBN, the contribution of 
GWs is characterized by as follows; 


Prad(T) -b pgw('T) 


30 

which leads to 


2 + — {4 -b 2(A^i, -f AAkjw)} 


(7) 


Pgw(T) 

= Prad(T) X ^ X 2 X A7 Vgw(T)/ |2 + ^(4 + 2Af,)| 

- ^Prad(T)AiVGw(T). (8) 

Noting 

GGw(Ti„) - Vh, (9) 


AA^upper- As is already mentioned, here it is assumed 
that any physical mechanisms, both known and un¬ 
known, contribute positively to A^eff- However, it would 
be worthwhile to note that at least there are examples 
where A"eff decreases (lol - [r^ . With this in mind, the re¬ 
quirement above, ANqw < AiVupper, is translated into 
an upper bound on Vh from (flOll as follows: 

,n, 

On the other hand, at the epoch of photon decoupling. 


Prad(T) -b PGw('T) 
^2 


TT‘ 


{2 + 2 X ^ ^ iN. + AiVGw)| , (12) 


which yields 

2x7 

11gw(T) =- 1 A!-x -AA^gw — O.ISAA^gw- 


2 + 2x|(A)^/3iV, 
So in this case we find 


Vh < O.I3AiV, 


upper 


/ PS (Tin) 
V Ps(T) 


1/3 


(13) 

(14) 


Note that this constraint depends on 35 , which one may 
regard as a drawback of these methods since gs is uncer¬ 
tain especially at high temperatures. It also depends on 
other potential entropy productions [i^. On the other 
hand, the PBH constraint does not much depend on gs 
nor other entropy productions. Assuming the standard 
model of particle physics, we use the following for sim¬ 
plicity: 


{ 100, (150MeV < T) 

gs{T) - 11, (IMeV < T< 150MeV) (15) 

[ 4 , (T < IMeV) 

Note that T ~ 150MeV corresponds to the QGD phase 
transition, and T ^ IMeV to the electron-positron anni¬ 
hilation. In order not to spoil BBN, we set, following (d^ . 
AVupper = 1-65 as a 95% C.L. upper limit, which can 
constrain GWs only of the scales smaller than the comov¬ 
ing horizon at the time of BBN, namely, 70kpc“^ < k. 
Noting the comoving scales which enter the horizon at 
T ~ 150MeV and T ^ IMeV correspond to fc ~ 4pc“^ 
and k ^ lkpc~^ respectively, we obtain the following 
constraints 


ANgw{T) can be written as 

AVgw(T) = |ggw(T) ^ yT. (^1^) . (10) 

In the literature an upper bound on AA"eff < AVupper 
is translated into an upper bound on AVgWi AA"gw < 


f 0.56, (4pc ^ < k) 

\ 0.27, (70kpc”^ <k< dpe"!) 


(16) 


As for CMB constraints, in Q the use of homogeneous 
initial conditions of GWs is advocated for SGWB gener¬ 
ated, for instance, by quantum fluctuations during infla¬ 
tion. In this case 95 % upper limits are AVupper = 0.18 
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(4^ , which correspond to 


r 0.068, 

(4pc ^ ^k) 


Vh < \ 0.033, 

(lkpc“^ ^ k< 4pc“^) 

(17) 

[ 0.023, 

(lOGpc"^ <k< Ikpc”^) 


For adiabatic initial conditions of GWs we refer to 


Aeff = 3.52lj];^| (95%; PZancfc-kWP-khighL-kHo+BAO) 

of [I| to set AiVupper 

= 1.00 [|: 

(18) 

r 0.38, 

(4pc“^ < k) 


Vh < < 0.18, 

(lkpc“^ ^ k < 4pc“^) 

(19) 

[ 0.13, 

(10Gpc“^ ^ k< lkpc“^) 



One may not regard these constraints as meaningful, be¬ 
cause upper limits correspond to the amplitude of GWs 
which is (almost) non-linear. 

The current energy density of SGWB, Oqw.Oj is 
also constrained by LIGO and Virgo, most severely in 
the band 41.5 — 169.25Hz as ^Jgw.o ^ 5.6 x 10“® x 
log(169.25/41.5) ~ 8 X IQ-^ [|. ’ Noting Hgw,o ~ 
z~q (4/100) ~ 10~'^Vh (-Zeq 3000 is the redshift 

at the matter-radiation equality, and the factor re¬ 
flects Ogw oc (l-|-z)/(l-|-Zen) during a matter-dominated 
era), we have Vh < 0.08 0 ■ 

Pulsar timing arrays (PTAs) have also been used 
to constrain GWs. Following [4^ we use the most 
stringent upper bound around / = 5.72xl0“®Hz 
4 X lO^Mpc-^), Qgw Ze-i(4/ll)i/37?^ < 2 x lO-®, 
which leads to < 1 x 10“^. 

Note that GW detectors or PTA experiments usually 
target GWs on a relatively limited frequency range, while 
cosmological methods like PBHs probe primordial GWs 
on a wide range of frequencies, and this is another ad¬ 
vantage of PBHs in constraining primordial GWs (see 
FigH]). 

Discussion.- As already mentioned PBH constraints 
are applicable from ^ Gpc all the way down to ^ 0.3m 
if we assume the number of e-foldings during inflation is 
sixty. Note that the exclusion of an overproduction of 
smallest PBHs (Mpbh ^ 10^(?) depends on the assump¬ 
tion that stable Planck mass relics are left over at the 
end of Hawking evap oration, which contribute to cold 
dark matter (see [i^, [l^ and references therein). The 
range of comoving scales corresponding to Mpbh 
is roughly < 50m, and so PBH upper bounds in this range 
depend on this assumption. If Planck mass relics are not 
formed, to what extent an overproduction of PBHs with 
MpBH lO^g is cosmologically problematic is uncertain. 
Such an overproduction of smallest PBHs may lead to an 
early matter-dominated era, during which PBH binaries 
are formed and emit GWs, or larger PBHs may form due 
to merger taking place after the collapse of perturbations 
of PBHs’ density, thereby leaving observable traces (i^ . 
Therefore, in principle one may still exclude such an over¬ 
production of smallest PBHs even without the left over 
of Planck mass relics to fully validate our upper bounds 


on smallest scales, though we do not discuss it in detail 
here. 

We can also constrain a blue tensor power spectrum of 
the following form; 

Vh(k) = rT’c(fcref) ( ^ j ’ < fcmax), (20) 

where is the dimensionless power spectrum of the 
curvature perturbation on uniform-density hypersurfaces 
and r is the so-called tensor-to-scalar ratio, here de¬ 
fined at fcref= O.OlMpc”^. Let us assume r = 0.2 and 
'P(kref) = 2.2 X 10“® following and also fcmax 
lGpc“^e®° ^ 3m“^. If we simply require Thikma,^) ^ 
0.04, we obtain ut < 0.3, which is tighter than other 
constraints such as BBN constraints shown in . Once 
more, this PBH bound does not depend on the details of 
a potential early matter-dominated era phase nor an en¬ 
tropy production, on which constraints other than PBHs 
are sensitive (4^ . If we use a more secure but weaker 
constraint on PBHs of around 10^5 by their entropy pro¬ 
duction [T^, we may require 'P^(20km~^) < 0.4, leading 
to nx ^ 0.4, still tighter than other constraints. To con¬ 
clude, PBHs can provide important constraints also on a 
blue spectrum, though a more careful analysis would be 
merited. 

We have imposed < 0.l5r,th to obtain PBH con¬ 
straints, eq. ©■ The factor of 0.1 here is to ensure the 
sufficient rareness of formation of PBHs. Note that to 
discuss the probability of PBH formation one needs to 
know the PDF of induced radiation density perturba¬ 
tion, which is in general highly non-Gaussian since it is 
generated by second-order tensor perturbations. In Q 
we assume Gaussianity of tensor perturbations on top of 
the delta-function power spectrum, and numerically con¬ 
firmed the factor of roughly 0.1 is indeed necessary to 
avoid an overproduction of PBHs (compare eq.© and 
the PBH constraints in FigH]). Strictly speaking, obser¬ 
vational constraints on PBHs depend on their mass [l^ 
and so how rare PBH formation has to be also depends 
on their mass, hence the scale dependence of PBH con¬ 
straints in FigH] HI . 

This argument implies PBH constraints on tensor per¬ 
turbations depend on the statistics of tensor pertur¬ 
bations, determining the statistics of induced density 
perturbation, just as PBH constraints on scalar per¬ 
turbations depend on the statistics of scalar perturba¬ 
tions [sO] . If high-cr realizations of tensor perturbations 
are suppressed (enhanced) in comparison to a Gaussian 
case, PBH constraints on tensor perturbations are tighter 
(weaker). 
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